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Super-Kamiokande has reported the results for the lepton events in the atmospheric neutrino 
experiment. These results have been presented for a 22.5kT water fiducial mass on an exposure 
of 1489 days, and the events are divided into sub-GeV, multi-GeV and PC events. We present 
a study of nuclear medium effects in the sub-GeV energy region of atmospheric neutrino events 
for the quasielastic scattering, incoherent and coherent pion production processes, as they give the 
most dominant contribution to the lepton events in this energy region. We have used the atmo- 
spheric neutrino flux given by Honda et al. These calculations have been done in the local density 
approximation. We take into account the effect of Pauli blocking, Fermi motion. Coulomb effect, 
renormalization of weak transition strengths in the nuclear medium in the case of the quasielastic 
reactions. The inelastic reactions leading to production of leptons along with pions is calculated in 
a A- dominance model by taking into account the renormalization of A properties in the nuclear 
medium and the final state interaction effects of the outgoing pions with the residual nucleus. We 
present the results for the lepton events obtained in our model with and without nuclear medium 
effects, and compare them with the Monte Carlo predictions used in the simulation and the experi- 
mentally observed events reported by the Super-Kamiokande collaboration. 

PACS numbers: 12.15.-y,13.15+g,13.60Rj,23.40.Bw,25.30Pt 



I. INTRODUCTION 



There are now many evidences that neutrinos oscillate and the neutrinos are not massless. These come from 
the experiments performed with atmospheric Ql-Qj reactor solar [lol l fill neutrinos and neutrinos obtained 

from the accelerators as in the experiments performed by MiniBooNE [I2||-jl6|. SciBooNE p7| . K2K [l^-[22| and 
MINOS collaboration [2^- [2^. The information from these experiments puts the limits on the solar and atmospheric 
neutrino mass differences viz. Aml^i^^ = 7.65toll x IQ-'^eV'^ and |AmatmP = 2.4lo;i[3 x lO'^eV'^. Two out 



of the four parameters of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix i.e. sm^{9i2) — 0.304lg'Q^g and 
sm^(2^?23) = 0-50lo Q6 known with a better accuracy while sin^{6i3) and CP violating phase 6 are still 

unknown. There is some information on 5171^(6*13) (< 0.01 at la), but there is no information on S. With these 
experimental observations lot of theoretical as well as experimental activities are going on. Now the aim of the 
experimentalists is to determine with better precision the various parameters of the PMNS matrix, absolute masses 
of the different flavors of neutrinos, to see whether the neutrino mass hierarchy is normal or inverted and the CP 
violation in the neutrino sector exists or not. On the theoretical side recently it has been emphasized in a series 
of neutrino workshops and conferences like Nuint [H], NuFact [H], NOW [Hi, etc. that at the neutrino energies 
of a few GeV energy region the study of neutrino nucleus cross section is very important which is relevant for the 
experiments with atmospheric neutrinos at Super-KamiokandeQ, ICARUS [soj . accelerator neutrinos at MiniBooNE, 
K2K, T2K [3]|- [3^, N OvA [s^ and the future experiments planned with Beta beams MhlHl and superbeams at 
Neutrino Factories [33|,[3l|. Neutrino-Nucleus cross section is one of the inputs in predicting the neutrino event rates. 
In the neutrino nuclear scattering process nuclear medium effects should be taken into account for writing the neutrino 
generator Monte Carlo codes which are used in analyzing the neutrino oscillation experiments. These oscillation 
experiments use various nuclear targets like ^^C, ^^O, ^'^Ar, ^^Fe, etc. For example, Super-Kamiokande(SuperK) is 
a 50kT water Cerenkov detector observing neutrinos from the terrestrial as well as accelerator neutrino sources. The 
nuclear target is ^^O in water(i?20), and neutrinos (antineutrinos) are interacting with the free protons as well as with 
the nucleons inside the oxygen nucleus. The Monte Carlo simulation of the lepton events uses Smith-Moniz model 
[3^ for the quasielastic process which does not include the effect of nuclear medium arising due to nucleon-nuclcon 
correlations but includes only the effect of Pauli principle and Fermi motion in a Fermi gas model. In the case of 
inelastic reaction for incoherent and coherent pion production it uses Rein and Sehgal model [40l . l4l| with the inclusion 
of the nuclear effects arising due to final state interactions of pions with the nucleus like pion absorption and pion 
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scattering. Furthermore, there are also some recent calculations on multipion production and deep inelastic neutrino 
reactions showing that the nuclear effects can be important in these energy regions. 

In this work, we have studied the lepton event rates for the atmospheric neutrinos at SuperK with and without 
the nuclear medium effects and compared our results with the experimental observed events and also with the Monte 
Carlo predictions for the events used by SuperK collaboration These results have been presented for a 22.5kT 
water fiducial mass on an exposure of 1489 live days, and we have taken the sub-GeV events in our analysis which have 
been classified as the region in which lepton's energy Ei < 1.33GeV and minimum observed momenta of electrons 
and muons are lOOMeV and 200MeV respectively [2}. Our Monte-Carlo analysis of the events has been done by 
considering the nuclear medium effects in the quasielastic, incoherent and coherent pion production processes, as they 
give the most dominant contribution in the sub-GeV energy region of atmospheric neutrino events. We have used the 
atmospheric neutrino flux given by Honda et al. 

In the case of quasielastic reaction, the effects of Pauli principle and Fermi motion are included through the Lindhard 
function calculated in a local density approximation. The renormalization of the weak transition strengths is calculated 
in the Random Phase Approximation(RPA) through the interaction of the p-h excitations as they propagate in the 
nuclear medium using a nucleon-nucleon potential described by pion and rho exchanges. The single pion production 
in the sub-GeV region is dominated by the resonance production in which a A resonance is excited and decays 
subsequently to a pion and a nucleon. When this process takes place inside the nucleus, there are two possibilities i.e. 
the target nucleus remains in the ground state leading to coherent production of pions or is excited and/or broken 
up leading to incoherent production of pions. We have considered both the production processes in the A resonance 
model in the local density approximation to calculate single pion production accompanied by a lepton from the oxygen 
nucleus. The effect of nuclear medium on the production of A is treated by including the modification of A properties 
in the medium. Once pions are produced, they undergo final state interactions with the residual nucleus. We have 
taken the final state interaction effects for both the incoherent and coherent pion production processes. This work 
is based on our study of the nuclear medium effects in the neutrino (antineutrino) induced reaction on the various 
nuclear targets like ^^C, ^^O, "^^Ar, ^^Fe etc., in the local density approximation, which have been applied to low and 
intermediate energy neutrinos for the charged current quasielastic process and the inelastic pion production process 
for the incoherent and coherent lepton production accompanied by a pion [43|-[50|. 

The plan of presentation is as follows. In section-2, we describe the charged current neutrino(antineutrino) in- 
duced quasielastic inclusive production of leptons from the nucleus. In section-3, we describe the charged current 
neutrino (antineutrino) induced inelastic production of leptons accompanied by a pion from the nucleus, where we 
describe the incoherent pion production as well as coherent pion production processes. In section-4, we present our 
results for the total scattering cross section as well as distribution averaged over the atmospheric neutrino flux 
given by Honda et al. [l^l,!^?! for the SuperK site. Furthermore, we discuss the dependence of the cross sections and 
distribution on the axial dipole mass Ma and the various parameterization of the form factors discussed in the 
literature recently for the isovector form factors in the case of quasielastic scattering and N-A transition form factors 
in the case of inelastic scattering. We have averaged the total scattering cross section a{E) and Q^-distribution over 
the atmospheric neutrino flux given by Honda et al. for the SuperK site to obtain the total lepton production 
event rate and the results have been compared with the observed numbers at SuperK and also with the numbers used 
by them in their Monte Carlo 0. In section-5, we conclude our findings. 



II. QUASIELASTIC REACTION 



The basic reaction for the quasielastic process is a neutrino interacting with a neutron inside the nucleus is given 

by 

ly,{k) + n{p)^l-{k')+p{p');l = e-,^l- (1) 
The invariant matrix element for the charged current reaction of neutrino, given by Eq.([T]) is written as 



M = ^cosec 
v2 



(2) 



where Gp is the Fermi coupling constant (=1.16639x10 ^GeV dc{= 13.1°) is the Cabibbo angle and the leptonic 
weak current is given by 

- S(fc')7M(l-75)w(fc) (3) 

is the hadronic current given by 



u{p) 



(4) 
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where, ~ [k ~ k'Y is the momentum transfer square and M is the nucleon mass. F^2{'1^) the isovector vector 
form factors and F^(q^), Fp{q^) are respectively the isovector axial vector and pseudoscalar form factors. 

Using the leptonic and hadronic currents given in Eq. ([3]) and Eq. ^ , the matrix element square is obtained by using 
Eq.©: 



\M\' 



G 



f cos' 9c L(;Jj^-^ 



L^ii'J is the leptonic tensor calculated to be 



= EE;//, = L^,('')+zL;t.^-), where 
L^^.u^"^ = 8 [k^X + k'^k^- g^^ k ■ k'] and 



(5) 

(6) 
(7) 



For antineutrino induced reaction vi^k) + p{p) + n{p');l = 6"*",/!+, L^/^J is given by 



The hadronic tensor J'^" is given by: 



1 



-Tr 



(i)' + M)p'(^ + M)r'' 



(8) 



where 



{q')j^ + F^ {q^T^.j^ + FX (g')7M75 + F^ (9^)5^75 



(9) 



and f'' = 7O r-'t 7O 

The hadronic current contains two isovector vector form factors ^"1^2 (9^) '^^ the nucleons, which are given as 



Fr2iq')-Fl,iq')^Fl[,iq') 



(10) 



where F^^"'\q'') and F2''"' {q'^) are the Dirac and Pauli form factors of proton(neutron) which are in turn expressed in 
terms of the experimentally determined Sach's electric ^^"(q^) and magnetic G^j^^{q') form factors of the nucleons 
given by 



pP(")/'„2 



G^^{q') = Fr{q')-Fr{q') 
G^iq') = Fr{q') + ^Fr{q') 
This results in the following form of the isovector vector form factors ^"1^2(9^) to be used in Eq.(l4|) 



Friq') = (1- 



4M^ 



GTiq') 



AM 



2 G^Piq") 



F^iq'' 



4) ' [GlP{q')-G^-{q')] 



AA'P / 



(11) 
(12) 

(13) 
(14) 



G^E^{q') and G^£P{q') are described by Galster parameterization [5l| with Q' = —q' 



G|(Q2) ^ Gd(Q'), G'kiQ') = -Tfi,,GD{Q^)^n, Gd{Q^) = 



l + Xr. 



A,, = 5.6, T 



4M2 



4M2 



9M9^ = GniQ'), ^^kiO^^GniQ') 



(15) 
(16) 
(17) 
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with proton and neutron magnetic moments as /ip=2.79285/XAr and /i„=-1.913/iAr, respectively, ~ —q^ and 
My=0.84GeV. 

The isovector axial form factor is parametrized as 



1 + ^ 

Ml 



(18) 



and is obtained from the quasielastic neutrino and antineutrino scattering as well as from pion electroproduction 
data. We have used axial charge i^_A(0)=-1.267 and the axial dipole mass M^=l.lGeV, which is presently being used 
in the SuperK analysis [s^l- The world average value for M^=1.026±0.020 GeV [5^, which is consistent with the 
recent NOMAD results for Myi=1.05±0.02±0.06 GeV [H^l obtained from the quasielastic and reactions with 
carbon, however, the values for axial dipole mass obtained by the experiments performed at K2K and MiniBooNE 
differ from this value. The values reported by the K2K experiments are M^=1.14±0.11GeV[T^ from the SciBar 
detector and A/A=l-20±0.12GeV [l^ by the SciFi detector experiments and M4=1.23±0.2GeV recently reported by 
the MiniBooNE collaboration [l3|. 

The pseudoscalar form factor (Q^) is dominated by the pion pole and is given in terms of F^{Q^) using the 
Goldberger-Treiman relation as 

Recently several new parameterizations for electromagnetic isovector form factors (55j-[59| have been presented which 
are obtained from the fits to the electron scattering data. To see the dependence of the cross section on the various 
parameterizations of the electromagnetic form factors we have used the parameterizations given by Budd et al. (ssj 
known as BBA-03, Bradford et al. ^] known as BBBA-05, Bosted [HH as well as the parameterization given by 
Alberico et al. [1^ 

The form of the electric and magnetic Sach's form factor given by Budd et al. ^] (BBA03) for the nucleon is 



fP (q'^\ — }: (20) 

J _i_ 3 253Q2 _,_ 1 422Q4 ^_ 0.08582Q6 + 0.3318Q8 - 0.09371Q1" + 0.01076gi2 ^ ' 



Hp 1 + 3.104Q2 + 1.428Q4 + 0.1112Q6 - 0.00698ig8 + 0.0003705Q1" - 0.7063£;-05Qi2 

m( 



GUQ^) 1 



fin 1 + 3.043Q2 + 0.8548Q4 + 0.6806Q6 _ 0.1287g8 + 0.008912Qio 

The form of electric and magnetic Sach's form factor given by Bradford et al. [s^ (BBBA-05) is 

^eW ) - 1 + 11.1^+ 13.6^2 + 33.0r3 
GljiQ^) _ 1 + 0.150T 



Hp l + ll.lr+19.6r2 + 7.54r3 

2. 1.25r+1.30r2 



QTi cn^) = 

' 1 - 9.86r + 305r2 - 758r3 + 802t4 

GliiQ^) ^ l + 1.81r 

Hn ~ l + 14.1r + 20.7T2 + 68.7T3 

The parametrization given by Bosted [11] is 

^^^^^ ^ " 1 + 0.62Q + 0.68Q2 + 2.80Q3 + 0.83g4 ^^^^ 

GIAQ') 1 



Hp 1 + 0.35Q + 2.44Q2 + 0.50Q3 + 1.04Q4 + 0.34Q5 

m( 



GliiQ') 1 



Hn 1-1.74Q + 9.29Q2_7.63Q3 + 4.63Q4 

-1.25^nTGzj(Q2) 

1 + 18.3t 



GliQ') = 
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Alberico et al. ^59|] have parameterized recently the electromagnetic form factors of the nuclcon based on the recent 
experiments performed at Bates, MAMI and JLab. The form of electric and magnetic Sach's form factor given by 
them (5§| is 

(23) 







1 - 0.14t 




1-1 


- 11.18T+ 15.18t2 


+ 23.57r3 






1 + 1.07t 






1-1 


-12.30r-h 25.43r2 


+ 30.39r3 






l-|-2.13r 






1-1 


-14.53r + 22.76r2 


4- 78.29r3 


GUQ^) 




-0.10 


0.10 


(1- 


f2.83g2)2 1 (1 + 


0.43Q2)2 



Using these parameterizations of the isovector form factors discussed above, we calculate the hadronic tensor given by 
Eq.® and the matrix element square using Eq.([5]). With this |A^|2 the charged current quasielastic lepton production 
cross section is calculated. 

The general expression for the differential cross section for the reaction shown in Eq.([T]) is given by 

f Glcos'ec \ i2nyS^{k+p-p'-k') d^k' d^p' ^^^^ 
= [-^—) 4V(fc ■ k'r - mlAQ [2.nK, (2.)32A, LK^^/)^^-^^ (24) 

where f is the number of fermions in the final state. The double differential cross section cro(i?ej l^'l) for the basic 
reaction is then written as 

^ 2 

a,{Eu \k'\) = ^l;f tnM\H[q, + E^- E,] (25) 

In a nucleus, the neutrino scatters from a neutron moving in the finite nucleus of neutron density Pni^), with a local 
occupation number n„(p, r). In the local density approximation the scattering cross section is written as 

a{Ei,\k'\) = j2dvdp^^nn{p,r)aoiEi,k') (26) 

where cro{Ei, \k'\) is given by Ea. (j25p . The neutron energy En and proton energy Ep are replaced by -Endp]) and 
Epijp + q\), where p is now the momentum of the target neutron inside the nucleus. Inside the nucleus the neutrons 
and protons are not free and their momenta are constrained to satisfy the Pauli principle, i.e., p„ < pp^ and p'p{= 

\Pn + q|) > PFp, where pf„ and ppp are the local Fermi momenta of neutrons and protons at the interaction point in 

i_ 1 

the nucleus and are given by pPr^ — [37r2p„(r)] and pp^ — [3vr2pp(r)] ^ , pn{r) and Pp{r) are the neutron and proton 
nuclear densities which are given in terms of the nuclear density of the oxygen nucleus 

Pn{r)^^-^-j^p{r); pp{r) ^ jp{r) (27) 

where p(r) is the density of oxygen nucleus taken to be 3 parameter Fermi(3pF) density and the density parameters 
have been taken from Ref. ^6^. There are also other parametrizations for the nuclear (^^O) density available in the 
literature like Harmonic Oscillator (HO) density ^3] and Modified Harmonic Oscillator (MHO) [io] density and we 
have also studied the dependence of the cross section on the various nuclear densities. Furthermore, in nuclei the 
threshold value of the reaction i.e. the Q-value of the reaction((5r) has to be taken into account, which we have taken 
to be the value corresponding to the lowest allowed Fermi transition. 

These considerations lead to a modification in the 6 function used in Ea. ((25|) i.e. 6[qo + En — Ep] is modified to 
6[qo + En{p) — Ep{p + q) — Qr] and the factor 

j^nn{p, r)-^<5[go + En - Ep\ (28) 

occurring in Ea. (|26p is replaced by — (l/7r)ImC/jv('?0: 9)7 where UNiqo,<l) is the Lindhard function corresponding to 
the particle hole(ph) excitation|43| shown in Fig.(IT]) and is given by 

dp MnMp Unjp) [1 - np{p + q)] 
(27r)3 EnEp qo + En{p) - Ep{p + q) + ie 

where qQ=E^^ — Ei — Qr- For the antincutrino reaction the suffix n and p will get interchanged 
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FIG. 1: Diagrammatic representation of the neutrino self-energy diagram corresponding to the ph-excitation leading to vi+n — > 
+p in nuclei. In the large mass limit of the IVB(i.e.Mvi' — > oo) the diagram 1(a) is reduced to 1(b) which is used to calculate 
\M\^ in Eq.©. 



The imaginary part of the Lindhard function is obtained to be [42 



ImUNiqo, q) = -^^^ [Ef, - A] (30) 
27r \q\ 



with g2 < 0, Ep^ - qo < Ep^ and — ^ — ^ < Ep^, where Ep^ = y Pf„'^ + , Ep^ ^ ^ pp^'^ + Mp^ and 

A = Max 



Ai„, Ep^ - go, i — 



With inclusion of these nuclear effects the cross section a{Ei,) is written as 

(j{E,) = -2GF^cos'e, [ "^^\^dr [ "^^^ k'dk' [^"^^'^ dQ^-^L^.J'^'' ImUN[E,, - Ei - Qr,q]. (31) 

Jr^irr Jk'^,„ "'Q™„. E^^El 

The outgoing lepton when comes out of the nucleus, its energy and momentum are modified due to the Coulomb 
interaction. The Coulomb distortion effect on the outgoing lepton has been taken into account in an effective momen- 
tum approximation(MEMA) [4 7], [61] in which the lepton momentum and energy are modified. In the local density 
approximation, the effective energy of the lepton in the Coulomb field of the final nucleus is given by: 

E,ff ^Ei + T4(r), 

where 

V,[r) = Zf-aiTT (- r ^I^r'^dr' + H ^I^r'dA (32) 
V Jo Zf J,, Zf J 

This leads to a change in the Imaginary part of the Lindhard function occurring in Eq. (|3ip 

ImUNiE,, -Ei-Qr,q)^ ImUNiE,, -Ei-Q^- V,{r), q) 

In the nucleus the strength of the electroweak coupling may change from their free nucleon values due to the presence 
of strongly interacting nucleons. Conservation of Vector Current (CVC) forbids any change in the charge coupling 
while magnetic and axial vector couplings are likely to change from their free nucleon values. These changes are 
calculated by considering the interaction of ph excitations in the nuclear medium in Random Phase Approximation 
(RPA) as shown in FigO The diagram shown in Fig[2]simulatcs the effects of the strongly interacting nuclear medium 
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FIG. 2: Many body Feynman diagrams (drawn in the limit Mw —» oo) accounting for the medium polarization effects 
contributing to the process «/; + n ^ + p transitions 



at the weak vertex. The ph-ph interaction is shown by the wavy line in Fig[2]and is described by the tt and p exchanges 
modulated by the effect of short range correlations. 

The weak nucleon current described by Eq.Q gives, in the non-relativistic limit, terms like FAfTT^ and iF2^2W'''+ 
which generate spin-isospin transitions in nuclei. While the term ii^2 3Xf''+ couples to the transverse excitations, 
the term Fa(tt^ couples to the transverse as well as longitudinal channels. These channels produce different RPA 
responses in the longitudinal and transverse channels when the diagrams of FigOare summed over. This is illustrated 
by considering the contribution of a term like Fa<t^ in Eq.Q. The leading order contribution of this term to the 
hadronic tensor J'-' in the medium is proportional to F\5ijImUN which is now split between the longitudinal and 
transverse components as 



F^Si-jlmllN 



Pa [qi^j + {^tj - qiqj)] i^^u, 



N 



(33) 



The RPA response of this term after summing the higher order diagrams like Fig[2]is modified and is given by J^^p^ 



1 - UnVi 1 - UNVt 



(34) 



where Vi and Vt are the longitudinal and transverse parts of the nucleon-nucleon potential calculated with tt and p 
exchanges and are given by 



Vl{q) 



P 



Vt{q) = ^ 
ml 



^2 



2 - o2 
p q 



(35) 



A^r = l.iGeV ^ Cp — 2, Kp = 2.5GeV, and rup are the pion and rho meson masses, and g' is the Landau-Migdal 
parameter taken to be 0.7 which has been used quite successfully to explain many electromagnetic and weak processes 
in nuclei [g^jjg^. Recently, in a work by Nieves et al. [13], g' has been taken as 0.63. We have studied the dependence 
of cross section on the Landau-Migdal parameter by varying g' by 10% i.e. taking g' — 0.63 and 0.77. 

This modified tensor J^^p^ when contracted with the leptonic tensor gives the contribution of the F\ term 
to the RPA response. The effect of the A degrees of freedom in the nuclear medium is included in the calculation 
of the RPA response by considering the effect of ph-Ah and Ah- Ah excitations as shown in FiglSJ^b). This is done 
by replacing [/at by Un — Un + Ua, where [/a is the Lindhard function for Ah excitation in the medium and the 
expressions for Un and Ua are taken from the Ref. [6^. The different couplings of N and A are incorporated in Un 
and Ua and then the same interaction strengths Vi and Vt are used to calculate the RPA response. These effects have 
been recently discussed by Nieves et al. [64| . 
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Thus, in the presence of nuclear medium effects, the total cross section a{E^), is written as 

a{E,) = -2Gf' cos2 9c [ r^dr [ k'dk' /"^""^ dQ^-J—L^.J^-'p^mUNlE,, - Ei - - K(r), -?] (36) 

"'r„,„ Jk'^i^ "'Q?„,„ ^il^ 

where J'^pj^ is the modified hadronic tensor when RPA effects are incorporated. 

III. INELASTIC PION PRODUCTION 

In the case of inelastic production process of leptons, the leptons are produced along with one or multiple pions. 
Around the energy region of 1 GeV, in the energy spectrum of atmospheric neutrino, the inelastic lepton production 
is dominated by the processes where lepton is accompanied by a single pion. These pion production processes take 
place mainly through the excitation of A which subsequently decay into a pion and a nucleon. When these processes 
take place inside the nucleus then either the target nucleus remains in the ground state where it does not change its 
identity leading to coherent production of pions or the target nucleus is excited and /or broken up leading to incoherent 
production of pions. Recently, two pion production and quasielastic hyperons [6^.[67j have been studied and are found 
to make small contribution to the lepton production in the energy region of present interest. 

The basic reaction for the charged current neutrino(antineutrino) induced one pion production (CCItt) in nuclei, 
is that a neutrino(antineutrino) interacts with a nucleon N. The various possible channels contributing for one pion 
production processes are 

Mk)+Pip) ^ /-(fc') + A++(P) (37) 

\ p + 7r+ 

iyi{k)+n{p) ^ r(fc') + A+(P) (38) 

\ n + 7r+ 

i?i{k)+p{p) ^ l+{k') + A°{P) (39) 

\p + 7r~ 
\n + 7T° 

Pi{k)+n{p) ^ l+{k') + A-{P) (40) 

\ n + 7r~ 

Now we shall present the formalism for the incoherent and coherent pion productions in brief. 



A. INCOHERENT PION PRODUCTION 



In the case of incoherent pion production in A dominance model, the weak hadronic currents interacting with the 
nucleons in the nuclear medium excite a A resonance which decays into pions and nucleons. The pions interact with 
the nucleus inside the nuclear medium before coming out. The final state interaction of pions due to elastic, charge 
exchange scattering and the absorption of pions leads to reduction of pion yield. The nuclear medium effects on A 
properties lead to modifications in its mass and width which have been taken from the work of Oset et al. (68j . 

In the case of incoherent one pion production process, the hadronic current J'' for the A excitation from the proton 
target is given by 

= V3*„(pa)0^'^«(p) (41) 

where ipaiPA) and u(p) are the Rarita Schwinger and Dirac spinors for the A and the nucleon, of momenta pA and 
p, respectively. O^^ is the N-A transition operator given by O^" = O^" + O^/ where 
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and 

- ^^{9''^A-q''r) + ^^ia''^1-PA-q>'^) + Ci{q')9-^ + ^^q^q'^ (43) 

where = P + 9; C'j^(i=3-6) are the vector and (7/^(1=3-6) are the axial vector transition form factors. 

The conserved vector current (CVC) hypothesis imphes Cg'(g^)=0. The other form factors CY {i = 3 — 5) are related 
in terms of the isovector electromagnetic form factors of the p — > A+ electromagnetic transition, and are determined 
from the analysis of data on photoproduction and electroproduction of A. 

The N-A vector transition form factors given by Schreiner and von Hippel[69l| are 

C^iQ^) = 2.05 (1 



V Ml 

CYiQ') - (44) 



with Ma as the invariant mass of the ttN system. 

The N-A axial vector transition form factors are given bv[6E 



C,^(Q^) = C^(0) (1 + ^) 7l--^^), ^ = 3,4,5 (45) 



-2 

Ml) V' W+W) 

with C^(0) = 0, C^(0) = -0.3, C^(0) = 1.2, ag ^ 63 = 0, 04 === 05 ^ -1.21, 64 ^ h =2 GeV^, = 1.05 GeV. 
While the various other parametrization for the N-A transition form factor has been discussed in literature 

The parametrization given by Lalakulich et al. [72] for the N-A transition form factors are given by 

Cr{Q')^Cr{0) (1 + -^) V, , z = 3,4,5. (46) 



where 



= 1 H TT for i = 3,4 and 

' AM?, ' 



-1 



V, = 1 + ^ for i = 5. (47) 



and 



QW) = C.» (l + f )" (l + l^)^ (48) 
Leitner et al. (t^I and Paschos et al. [zH use the following form of the N-A vector transition form factor 

where C|'(0) = 1.95, Cl'(O) = -# C^iQ), {Q^) = with W as the center of mass energy (V(P + a)^) and Ma 
as the mass of A. For the axial vector part they use the same parametrization as used by Lalakulich et al.fT^]. In 
the reactions given by Ea l571HDl a A is produced, which subsequently decays into a nucleon and a pion, for example 
interaction of a neutrino with a proton inside the nucleus is given by z^;(fc) -l-p(p) — > ' (fc') +p(p') + 7r+(fc7r), for such 
a process the transition matrix element M. fi is given by 

Mf, = V3^^^^u{p')KV.,0^'^IMp) (50) 
V 2 ^TL■^ 



10 



where is the leptonic current given by Eq.(I3]), O'^^ is the N-A transition operator given by Eq. (|l^ & Ea. (|15)) . and 
V^^ is the A propagator in momentum space which is given by : 



P2 _ Ml + iMaT ^^^^ 



with P"^ as the spin-3/2 projection operator given as 



sptns 

and the delta decay width F is taken from [g^I i.e.: 



|qcm| is the pion momentum in the rest frame of the resonance given by 

_ - ml - M2)2 - 4m2M2 

and M is the mass of nucleon and W is the center of mass energy. In the nuclear medium the properties of A hke its 
mass and decay width F to be used in Eg . ([STl ) are modified due to the nuclear effects. These are mainly due to the 
following processes. 

(i) In the nuclear medium As decay mainly through the A — > Nn channel. The final nucleons have to be above 
the Fermi momentum fc^? of the nucleon in the nucleus thus inhibiting the decay as compared to the free decay of 
the A described by F in Eq.((53|. This leads to a modification in the decay width of delta. We have taken these 
modifications given by |68] where the modified delta decay width F is given as 

t^TxFikp,EA,kA) (54) 
Here F{kF, E^, ^a) is the Pauli correction factor given by [68| : 

kA\(\cm\+EAE' ^EpW 



Ef = \/ M'^ + kp, fcA is the A momentum and Ea = \/W + k\. 

(ii) In the nuclear medium there are additional decay channels open due to two and three body absorption processes 
like AA^ — > NN and ANN — > NNN through which A disappears in the nuclear medium without producing a pion, 
while a two body A absorption process like AA^ — s- t:NN gives rise to some more pions. These nuclear medium effects 
on the A propagation are included by describing the mass and the decay width in terms of the self energy of A. These 
considerations lead to the following modifications in the width F and mass Ma of the A resonance. 

f f 

- ^ - - /ttiEa and Ma Ma + ReT,A- (56) 

The expressions for the real and the imaginary parts of Sa are [6^ : 

ReJ:A = AO—MeV and 
Po 

= C«(ii)%G„(ii)%C„(ii)' (57) 

In the above equation Cq accounts for the AA^ ttNN process, Ca2 for the two-body absorption process AA^ — > A^A^ 
and Ca3 for the three-body absorption process ANN — * NNN. The coefficients Cq, Ca2, Ca3 and a, (3 and 7 are 
taken from Ref. I68l. 
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Thus, in the local density approximation the expression for the total cross section for the neutrino induced charged 
current Ivr^ production from proton target is written as 

I rr-ma^ rQ'ina^ r^'ma^ + 1 

aA{E) = 7- — TT- / pp{r)dr I dQ^ x / dk' / dcosO-^ 
[^^r Jr„,^„ JQl.i^ "'-1 

^\k'\\K\ 1 



X 



MElEi e,^^e^^^_Wlc^^(^ 



SS|M/zr (58) 



Similar expression is obtained for the cross section for neutrino induced charged current Itt^ production from neutron 
target with pp replaced by ip„. A factor of ^ comes with p„ due to the Clebsch Gordan coefficient occurring in the 
production of tt'*' from the neutron target [v^ + n ^ + A+, A+ — )■ n + tt"*") as compared to the 7r+ production 
from the proton target. For the antineutrino induced Itt^ production process the factor of ^ will come with pp i.e. 
7r~ would get produced dominantly from neutron target. 

The pions produced in these processes inside the nucleus may rescatter, produce more pions or may get absorbed 
while coming out from the final nucleus. These are treated using Monte Carlo simulations by generating a pion of 
given momentum and charge at a point r in the nucleus. Assuming the real part of the pion nuclear potential to 
be weak as compared with their kinetic energies, they are propagated following straight lin es till they are out of the 
nucleus. At the beginning, the pions are placed at a point (r = b,Zi„), where Zin = — -\/R^ — |bp, with b as the 
random impact parameter, obeying |b| < R. R is upper bound for the nuclear radius, which is chosen to be such 
that p{R) ~ lO^'^poi with po is the normal nuclear matter density. The pion is then made to move along z-direction 
in small steps until it comes out of the nucleus. We have taken the results of Vicente Vacas for the final state 
interaction of pions which has been discussed in Ref. [ill . 

B. COHERENT PION PRODUCTION 

The Vfj_ induced coherent one pion production on ^^O target is given by z^^ O ^ p^ +1^ O + 7r+ for which the 
cross section is given by Ea. (|55|) . However, the matrix element Aifi is now given by 

Mf^^ ^cosej^J^J^iq- k^) (59) 
v2 

where is the leptonic current given by Eq. ((3|) and is the hadronic current given by [50l | 

J = ^kllA V uMKaV'^Ox^Urip) (60) 

paX gjygjj |-,y Eq. ([5T|) . is given by Eqs. and and u(p) is Dirac spinor for the nucleons. 
Here J-{q— fcvr) is the nuclear form factor, given by 



^(•f- K) = / d^r 



When pion absorption effects are taken into account using the Eikonal approximation then the nuclear form factor 
T{q — kyr) is modified to T{q~ k^), which is calculated in Eikonal approximation to be [75| : 



oc 







T{q-k^) = 2TT I bdbj dz p{h,z) Ja{klb) e'^^^-^^^-e-'f^''^''-'' (62) 

where 



f{b,z)= mpib,z'))dz' 
Jz 2\k^\ 

k\ and k\ are the longitudinal and transverse components of the pion momentum and H is the self-energy of pion, 
the expression for which is taken from Ref. [75| and is given by 

nwS.O) ^ \ .(MO ^^^^ («3, 

Using the matrix element given by Eq. (j59p and the modifications in the A mass and width in the nuclear medium 
given by Eq. ([55)) we calculate the total scattering cross section a given in Eq.(|5^ 
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FIG. 3: Total scattering cross section((T) as a function of neutrino energy for the charged current (a) Quasielastic (b) Incoherent 
and (c) Coherent processes for induced reaction in ^^O. 




FIG. 4: Total scattering cross section((7) as a function of neutrino energy for the charged current (a) Quasielastic and (b) 
Incoherent processes for 1/^, v^^, Ve and Ufj, induced reaction in ^^O. These results are presented for a calculated with RPA in 
CCQE scattering (Fig. 4(a)) and with medium and pion absorption effects in the case of CClvr production processes (Fig. 4(b)). 



IV. RESULTS AND DISCUSSION 
A. Total Scattering Cross Section a 

For the charged current quasielastic(CCQE) reaction, the numerical results are obtained from Eg. ([36]) using the 
expression for the form factors given by Bradford et al. [s^ with vector dipole mass My=0.84GeV and axial dipole 
mass Myi=l.lGeV. In the case of charged current induced incoherent and coherent pion productions the results for 
the total cross sections are obtained from Eq. ([58|l using the matrix elements given in Eq. ([50| & Eq. ([59| respectively 
and the N-A transition form factors given by Lalakulich et al. [t^] given by Eq. (j46|) and Eq. ((48|) with M^=l.l GeV. 
In FigsO^-c, we have shown the results with nuclear medium effects for the total cross section cr for CCQE and 
incoherent & coherent CCl7r+ production cross sections in "'^^O, for induced reaction. 

For the quasielastic process, in the case of charged current induced lepton production cross section the results 
have been presented for the cross section calculated for the free case, with nuclear medium effects without RPA i.e. 
our local Fermi gas model, and with nuclear medium effects including RPA. These results have been shown in Fig. ([3^). 
We find that when the cross section is calculated in the local Fermi gas model the reduction in the cross section is 
around 18% at E^^=0.4GeV and around 10% at Ei/^_=l-3GeV from the cross sections calculated for the free case. 
However, when we encorporate the RPA effects, there is further reduction in the cross section which is about 30% at 
E^^=0.4GeV, 15% at E^^^l.OGeY and around 12% at i;^^=2-3GeV. 

The numerical results for the induced incoherent Itt"*" production process have been shown in Fig.(I5j3) and we 
find that the nuclear medium effects lead to a reduction of around 40% for Ey^=0.4 GeV, 30% for Ei,^^=1.5 GeV and 
28% for E^^=3 GeV. When pion absorption effects are also taken into account along with the nuclear medium effects 
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FIG. 5: Total scattering cross section((T) as a function of neutrino energy for (i/^) induced reaction in ^^O. These results 
have been presented for the cross sections calculated by using different models. 
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FIG. 6: Total scattering cross section((T) as a function of neutrino energy for the charged current quasielastic lepton production 
process induced by (a) z/^ and (b) in ^^O in the local FGM with RPA effects. The present curves show the Ma dependence 
on the cross section. 



there is a further reduction in the cross section which is around 15% for Ey^=1.0 GeV and 12% for Ey^=3 GeV. 
For induced coherent Itt'^ production process the numerical results have been shown in Fig.([3J;) and we find that 
the nuclear medium effects lead to a reduction of around 40% for E^^^=0.8 GeV, 24% for Ey^_=2.0 GeV and 15% at 
Ey^=3 GeV. When pion absorption effects are also taken into account along with the nuclear medium effects there 
is a further reduction in the cross section which is around 50% for Ei,_=0.8 GeV, 35% for Ei,_=2.0 GeV and 25% 
E,^,=3 GeV. 

Thus, in the case of incoherent production of pions, the reduction due to nuclear medium effects in the production 
process is larger than the reduction due to final state interaction while in the case of coherent pion production, 
the reduction due to final state interaction is quite large as compared to the reduction due to the nuclear medium 
effects. Furthermore, the contribution of the cross section calculated in the case of coherent pion production process 
in the nuclear medium with final state interaction effects is around 6-7% to the total (incoherent -|-coherent) one 
pion production cross section in the energy region of OAGeV < E < 3GeV, and due to this we have not discussed 
the form factor dependence, distribution, etc. in the results presented here for the coherent pion production 
process as their contribution to the total lepton events in the sub-GeV energy region of present interest is not very 
significant. Our results for the coherent process also agrees with the other recent calculations performed by the various 
groups [2|,[73,[zi. 

In Figs.(|l^) and we have presented the results for the total scattering cross section cr as a function of neutrino 
energy for the charged current lepton production process induced by neutrino(antineutrino) in ^^O in the case of 
quasielastic and incoherent Itt production processes. These results have been presented for I'e, i^^, and i^^ in the 
local Fermi gas model with RPA effects in the case of quasielastic process and with nuclear medium and final state 
interaction effects in the case of incoherent Itt"*" (neutrino) and l7r~(antineutrino) production processes. 

In Figs. ([5^) and (03), we have compared our results for the and i?^ induced charged current quasielastic lepton 
production cross sections in ^^O, obtained in the local Fermi gas model with and without the RPA effects, with the 
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FIG. 7: Total scattering cross section((T) as a function of neutrino(antineutrino) energy for the charged current incoherent 
l7r"'"(l7r~) production from ^^O target with nuclear medium and final state interaction effects. The present curves show the 
Ma dependence on the cross section. 



results obtained in the Fermi gas model given by Smith and Moniz [39| and Llewellyn Smith [79|, which have been 
used in some of the Monte Carlo generators. We find that our results in the local Fermi gas model are in fairly good 
agreement (within 2%) with their results [s^, [t^ in the case of Vy, induced process, while in the case of process the 
results obtained in our local Fermi gas model are within 3-4% with the results obtained by using Llewellyn Smith's [7§| 
Fermi gas model, however, the results obtained by Smith and Moniz [s^ Fermi gas model is about 5-6% higher. We 
have also compared our results in the local FGM with the non-relativistic FGM of Gaisser and O' Connell [8^ (not 
shown here) and found the cross sections to be within 2-3 % at the neutrino energies of the present interest. When 
RPA correlation effects are taken into account the cross section decreases. We find that our results for the total 
scattering cross section uiE) in the case of charged current neutrino induced process calculated in the local FGM 
with RPA effects agree with the recent calculations performed by Leitner et al. [sij, Benhar et al.|82] and Nieves et 
al. H. 

We have shown in Figs.([n^) and ([BId), the effect of varying the axial dipole mass Ma on the total scattering cross 
section a for the and v^^ induced charged current quasielastic reaction cross sections in ^^O, obtained in the local 
Fermi gas model with RPA effects using Bradford et al. ^] parametrization of the weak form factors given in Eq. ((2T|) 
with vector dipole mass My=0.84GeV. We find that a increases by about 5% at £'i,^=0.5GeV and by 10% at E^^=2- 
SGeV when Ma is taken as 1.21GeV, while it decreases by about 3% at £'^^,=0.5GeV and by 5% at E^^=2-3GeV 
when M4=1.05GeV as compared to the cross sections calculated with Myi=l.lGeV. While in the case of induced 
reaction cross section the dependence on Ma becomes small which is around 6% at £'p^^ = l-3GeV when Ma is taken 
as 1.21GeV, while the decrease is about 3% at £'p^=l-3GeV when Af^=1.05GeV as compared to the cross sections 
calculated with M^=l.lGeV. 

We have studied the Landau-Migdal parameter {g') dependence(not shown here) given in Eq[35l expression for 
which has been used in calculating the cross section with RPA effects, on the total cross section. We find that a 
10% uncertainty in g' leads to a 5-6% of uncertainty in the cross section. Also we have studied the nuclear density 
dependence on the total scattering cross section(not shown here). Using other densities (modified harmonic oscillator 
or 3 parameter Fermi density) leads to 2-3% of uncertainty in the cross section. 

In Figs. ([7^) and ([7}d), we have shown the effect of varying the axial dipole mass Ma on the total scattering cross 
section a for the z^^ and induced Itt production cross sections in ^^O with nuclear medium and final state interaction 
effects. We have used the N-A transition form factor parameterizations given by Lalakulich et al. [T^I- We find that a 
increases by about 5% at E^=0.5GeV and by 12% at £^i,=2-3GeV when Ma is taken as 1.21GeV, while it decreases by 
about 3% at E,y=0.5GeV and by 6% at E,^=2-3GeV when MA=l-05GeV as compared to the cross sections calculated 
with Myi=l.lGeV. While in the case of i?^ induced reaction cross section the dependence on Ma becomes small which 
is around 10% at i<^p^=l-3GeV when Ma is taken as 1.21GeV, while the decrease is about 6% at Ei,—l-3GeY when 
MA=1.05GeV as compared to the cross sections calculated with M^=l.lGeV. 

Our results for the total scattering cross section a{E) in the case of charged current neutrino induced incoherent 
pion production process in the A dominance model calculated with nuclear medium and final state interaction effects 
agree with the numerical results of Leitner et al. 81] and Benhar et al.|82l]. 
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B. Differential Scattering Cross Section < > 



In this section we shall present the results for the flux averaged differential scattering cross section < > as 
function of Q^. This has 
al. |43j for the SuperK site. 



a function of Q^. This has been obtained by integrating over the atmospheric neutrino flux given by Honda et 



The flux averaged differential scattering cross section < > is defined as 

da rr m^mdE 



IE^^,„ dQ 

dQ2 ' (f>{E)dE 



< >= (64) 



where is the differential scattering cross section for the distribution and is the atmospheric neutrino 

flux. 

We discuss the nuclear medium modiflcation effects on Q^-distribution, the effect of Ma and different parameteri- 
zations of the various isovector form factors in the case of quasielastic process and N-A transition form factors in the 
case of Itt production process on distribution. In addition to these, we shall also present the results to show the 
dependence of distribution for the two fluxes Kaml997 and Kam2000 given by Honda [43j . 

In Figs.(IHti) and ((BJd), we present the results for the Q^-distribution in the case of charged current quasielastic 
lepton production process induced by electron and muon neutrino(antineutrino). The results have been presented 
for the Q^-distribution calculated in the local Fermi gas model with and without RPA effects. In the case of i^e 
and induced processes the results have been shown in Fig. ([8^). We find that in the case of I'e induced process 
the differential cross section is calculated in the local Fermi gas model with RPA effects the reduction in the cross 
section is around 42% in the peak region of (5^(=0.044GeV^) and around 30% at Q^=0.2GeV^ as compared to the 
cross section calculated without the RPA effects. For there is a shift in the peak region which is towards low 
(5^(=0.022GeV^) and the reduction is around 35% which is smaller than in the case of Ve induced process and at high 
(=0.1GeV^) the reduction is around 30% as compared to the cross section calculated without the RPA effects. 
When a cut on the electron's energy {Ee < 1.33G'ey) and momenta(pe > lOOMeV) is applied, then there is a small 
change in the spectrum which in turn leads to a small change in the event rates. 

In the case of induced process the reduction in the cross section when RPA effects are taken into account is 
around 40% in the peak region of Q^(=0.06GeV^) and around 30% at (5^=0. 2GeV^ as compared to the cross section 
calculated without the RPA effects in the local Fermi gas model. In the case of i/^, there is a shift in the peak region 
which is towards low (5^(=0.028GeV^) and the reduction is around 36% which is smaller than in the case of i/^ induced 
process and at high (=0.2GeV^) the reduction is around 22% as compared to the cross section calculated without 
the RPA effects. When a cut on the muon's energy {Ef^ < 1.33Gey) and momenta(pp > 2Q0MeV) are applied, then 
there is a large suppression in the distribution in the peak region as shown in the Fig|51 We find that the inclusion 
of RPA effects in our local Fermi gas model with a cut on muon's momenta and energy results in a large suppression 
in the event rates in comparison to the muon events calculated by using the local Fermi gas model without RPA 
effects and without applying cuts on the muon's energy and momenta. 

The distribution in the case of charged current electron and muon neutrino (antineutrino) induced incoherent 
l7r'^(l7r~) production processes have been shown in Figs. ((9^) and ([9}d). The results have been shown without and 
with the effect of nuclear medium as well as with the final state interaction of pions taken into account along with 
the nuclear medium effects. In the case of Ve and i^^ induced processes the results have been shown in Fig. ((2^). We 
find that in the case of Ve induced process when the cut is applied on the electron's energy (i^e < 1.33Gey) and 
momenta(pe > lOOAfe^), and the nuclear medium effects are also taken into account, the reduction in the cross 
section is around 36% in the peak region of Q^(—0.Q8GeV^) and around 33% at Q^=Q.7GeV^ as compared to the 
cross section calculated without taking nuclear medium effects into account. When final state interactions effects are 
also taken into account there is a further reduction of around 15% in the peak region and around 14% at Q^^O.TGeV^. 
For De, there is a shift in the peak region which is towards low Q^{—0.002GeV^) and the reduction is around 35% 
and at high (=0.2GeV^) the reduction is around 30% as compared to the cross section calculated without medium 
effects. Here we have also shown the results of the distribution calculated by including nuclear medium and final 
state interaction effects but without applying any cuts on lepton energy and momenta, it has been found that this 
results in an enhancement in the distribution particularly in the peak region of Q^. 

In the case of induced process, when the differential cross section is calculated by applying cuts on the lepton's 
energy and momenta, with the nuclear medium effects taken into account, the reduction in the cross section is around 
35% in the peak region of (5^(«0.1GeV^) as compared to the cross section calculated without the nuclear medium 
effects. When pion absorption effect is also taken into account there is further reduction of about 15%. In the case 
of Df^ induced process there is a shift in the peak region which is towards low Q^{~0.02GeV^) and the nature of 
reduction is almost the same as in the case of induced process. When there is no cut applied on the muon's energy 
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FIG. 8: < > vs for the quasielastic process induced by electron type(upper panels) and muon type (lower panels) (a) 



neutrino and (b) antineutrino in ^^'O. The results are presented for the Q^-distribution with and without applying cuts on the 
lepton's energy and momenta. The dashed(dotted) line is the result in the local FGM with(without) RPA effects and without 
cuts. The solid (dashed-dotted) line is the result in the local FGM with(without) RPA effects and with cuts {Ei < 1.33GeV^, 
Pe > lOOMeV and > 200MeV). 



and momenta, then there is a change in the nature of reduction. For example, in the case of induced process the 
reduction in the cross section calculated with the nuclear medium effects taken into account, is around 35% in the 
peak region of as compared to the differential cross section calculated without the nuclear medium effects, and 
when pion absorption effect is also taken into account there is further reduction of about 14%, while for i/^ induced 
process this reduction is around 30% in the peak region and with medium and final state interaction effects taken into 
account, the further reduction in the differential cross section is around 12%. To see the effect of applying cuts on 
the muon's energy and momenta, we calculate distribution with medium effects and final state interaction effects 
and find that when cuts are taken into account in the case of induced reaction the reduction is around 45% in the 
peak region and becomes 32% around =0.7GeV^ . However, in the case of induced reaction this reduction is 
around 35% in the peak region and becomes 24% around =0.2GeV^ . 

In Fig. pilk ). we show the results for the distribution in the case of charged current quasielastic lepton production 
process induced by z/g, P^, and in -'^^O calculated in the local Fermi gas model with RPA effects. The results 
for the distribution in the case of charged current incoherent l7r+ production process induced by i/g, v^, and Itt^ 
production process induced by Pg, in ^^O with nuclear medium and pion absorption effects have been shown in 
Fig.diniD). 

We have shown in Figs. pTk ) and pTb ). the effect of various parameterizations of the isovector vector form factors 
on the distribution in the case of charged current quasielastic lepton production process induced by and in 
^^O calculated in the local Fermi gas model with RPA effects. The results have been shown with the parameterizations 
given by Budd et al. [1^, Bradford et al. [56[, Bosted et al. ^58] and Alberico et al. [5§|. We find that the use of 
various parametrization for the isovector form factor results in a very small change in the distribution in the peak 
region. 

In Figs. p^ ) and (|12b ). the dependence of the various parameterizations of the N-A transition form factors on 
the Q^-distribution in the case of charged current l7r+ production process induced by z/g and Itt" production process 
induced by Vf, in ^^O with nuclear medium and pion absorption effects have been shown. These results are obtained 
by using the N-A transition form factors parameterizations given by Lalakulich et al. [T^I, Paschos et al. [7l| and 
Schreiner et al. [1^. We find that in the case of induced pion production process the differential cross section 
obtained by Paschos et al. I7l| is 5-7% smaller in the region of ~0.1-0.5GeV^, while the differential cross section 
obtained by Schreiner et al. yy\ is 5% smaller at low and which increases to around 10-16% for Q^=0.2-0.4GeV^, 
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interaction effects. The results for the distribution with nuclear medium & final state interaction effects but without putting 
cut on lepton's momenta and energy have also been presented here(dashed-dotted line). 
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FIG. 10: < > vs for (a) Quasielastic and (b)Incoherent processes induced by different fiavors of (a) neutrino and (b) 
antineutrino averaged over the atmospheric neutrino flux given by Honda et al. [4^.[43l]. The results have been presented after 
putting cut on lepton momenta and energy. 



than the cross section calculated with Lalakulich et al. f72l| parameterization. In the case of De induced process the 
use of various parametrization for the N-A transition form factors results in a very small change in the peak region 
of distribution. 

To show the dependence of the different fluxes at the Superkamiokande site for the Solar minimum and Solar 
maximum defined by Kaml997 and Kam2000 by Honda et al. [13], we have obtained the numerical results for the 
distribution in the case of charged current quasielastic lepton production process induced by and reactions in 
^^O calculated in the local Fermi gas model with RPA effects. The results for i^^ induced process is shown in Fig. p^ ) 
and for the i?^ induced process the results are shown in Fig. p3b ). We find that these two fluxes result in a very small 
difference in the distribution. Similarly in the case of inelastic one pion production cross section we find that the 
difference in the numerical results calculated using these two fluxes is very small. 
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Aluon Events 










Process 


Free Case 


r (oriVl 


rrr^T^/T "dtda 
r IjM W itn Kr A 


I 




3332 


2472 


1894 


11 




966 


620 


461 


III 


v^p ^^n[on free p due to H2) 


241 


24lf 


24lf 




P^ + 


4539 


3232 


2596 




Electron Events 










Process 


Free Case 


FGM 


FGM With RPA 


IV 


Ven e~p{m '-'^O) 


2332 


1754 


1278 


V 


UeP -> e+n(in '"O) 


609 


358 


266 


VI 


PeP — > e+n(on free p due to H2) 


152 


152t 


152t 




Ve + 


3093 


2264 


1696 



TABLE I: Total number of lepton events for a quasielastic process, f: For reaction on free protons the events would be the 
same in all the three columns. 





Muon Events 








Inelastic Process 


Free Case 


Medium effects 
with Pion abspn 


VII 


v^i,p — > ^^A^"^(on free p due to H2) 


154 


154 


VIII 


Ufj,p — > /i+A"(on free p due toH2) 


12 


12 


IX 


i^^'^0(/i~ accompanied by tv") 


171 


92 


X 


v^'O (/i* accompanied by vr" 


40 


23 


XI 


u^'^0{fi~ accompanied by 7r+) 


756 


409 


XII 


ul!^0{^'^ accompanied by tt^) 


179 


102 


XIII 


+ i^fi (Coherent) 


233 


30 


XIV 


+ (Quasielastic like events from Inelastic Process) 




344 






1545 


1166 




Electron Events 








Inelastic Process 


Free case 


Medium Effects 
with Pion Abspn 


XV 


VeP — » e~A'*'^(on free p due to H2) 


98 


98 


XVI 


UeP —* e+A"(on free p due to H2) 


6 


6 


XVII 


Ve^O (e~ accompanied by vr") 


99 


53 


XVIII 


i't'^O (e^ accompanied by tt") 


21 


12 


XIX 


ul^O (e^ accompanied by tt^) 


501 


269 


XX 


y^^O (e+ accompanied by ti~) 


91 


52 


XXI 


Ve + f'e (Coherent) 


148 


19 


XXII 


+ t^e (Quasielastic like events from inelastic process) 




200 




+ Ue 


964 


709 



TABLE II: Total number of lepton events for inelastic process. 



Process 


+ 




Free case(QE+Inelastic) 


4057 


6084 


FGM without RPA+Inelastic with nuclear medium 
and final state interaction effects 


2973 


4499 


FGM with RPA +Inelastic with nuclear medium 
and final state interaction effects 


2405 


3762 


Monte Carlo events 


2533.9 


3979.7 


Reported by experiments 


3353 


3227 



TABLE III: Total number of lepton events calculated in our model and its comparison with the observed lepton events by 
SuperK collaboration and the Monte carlo number used by them Qj. 
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FIG. 11: < ^^j- > vs for the quasielastic process induced by electron type (a) neutrino and (b) antineutrino with various 
parameterization of the isovector form factors, with Mv=0.84GeV and M^^l-lGeV. 




FIG. 12: < > vs for incoherent process induced by electron type (a) neutrino and (b) antineutrino averaged over the 
atmospheric neutrino flux with various form factors. 



C. Total lepton events 



Here we are going to present the results for the total number of lepton events for the sub-GeV energy region. These 
results have been presented for a 22.5kT water fiducial mass for 1489 days, and we have put a cut on the lepton's 
energy Ei < 1.33GeV and momenta of electrons and muons as Pe >100MeV and >200MeV. We have integrated 
the total scattering cross section a over the atmospheric neutrino flux given by Honda et al. [43| for the SuperK cite. 




FIG. 13: < ^^2- > vs for quasielastic process induced by muon type (a) neutrino and (b) antineutrino averaged over the 
atmospheric neutrino flux for the solar minimum and solar maximum given by Honda et al. [i^], [U- 
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The flux averaged cross section is defined as 

< (7 >= / a{E)(l){E)dE (65) 

where cr(E) is the total scattering cross section and 4'{E) is the atmospheric neutrino flux. 

In Table-I the lepton event rates have been obtained for the CCQE processes induced by i^i and Di {I = e,/i). In 
Table-II the lepton event rates have been obtained for the CCItt production due to i>i and i?i (l — e, /i)induced reactions 
in ^^O, as well as the leptons obtained from i^i and Di (l — e,/i)induced reactions on the free protons and the leptons 
accompanied by 7r° in the neutrino(antineutrino) induced processes. In Table-Ill, we have presented the total lepton 
events from vi + vi (e or /x) induced quasielastic and inelastic pion production processes. Here we have compared our 
final results (leptons obtained from CCQE reaction in local FGM with RPA effect + CCItt production with nuclear 
medium and final state interaction effects) with the experimentally observed lepton events by SuperK collaboration 
and also with the lepton events used in the Monte Carlo analysis of these events by the SuperK collaboration 2] . 

In our calculations for predicting the lepton events, we have considered the following channels. In the quasielastic 
process the contributions to the lepton events have been taken from the channels (i) i/;n — > Z~p(m ^^O), (ii) i'lP —> 
l~^n{in ^^O) and (iii) Dip l^n (on free p due to H2 )■ In the case of incoherent pion production process the various 
channels contributing to the lepton events are (i) vip /~A++ (on free p due to H2 ), (ii) vip Z+A^ (on free p 
due to H2 ), (iii) vi ^^O {l~ accompanied by tt" ), (iv) vi ^^O (/+ accompanied by 7r° ), (v) vi ^^O ( l~ accompanied 
by 7r+ ), (vi) vi ^^O ( accompanied by tt"). In the case of coherent process the contributions would come from 
(i) VI O r + 7r+ O and (ii) Di +^6 Q ^ 1+ + t^- O processes. 

In the case of incoherent and coherent pion production processes, when pion absorption effects are taken into 
account, the pions which are produced but get absorbed while coming out of the nucleus, the reaction produces only 
lepton. Such reactions give leptons without the pions in the final state and are labelled as quasielastic like events. We 
have also considered quasielastic like events in calculating the event rates. The effect of nuclear medium effects on 
the total number of lepton events can be summarised by saying that it leads to a reduction of 40% in the event rate. 



V. CONCLUSIONS 



In this work we have studied the nuclear effects in the charged lepton production in water induced by atmospheric 
neutrinos at the SuperK site. The energy dependence for the total cross section for the quasielastic, incoherent and 
coherent pion productions have been calculated in the local density approximation for neutrino (antineutrino) induced 
reactions in ^^O. We find that nuclear medium effects play a very important role in the study of a as well as fiux 
averaged differential cross section < > . The total lepton production event rate is compared with the experimental 
observed numbers at SuperK and also with the numbers used in their Monte Carlo[2]- 

We conclude the following: 

(i) In the case of neutrino induced charged current quasielastic lepton production process, the nuclear medium 
effects like Pauli blocking, Fermi motion effects, renormalization of weak transition strengths in the nuclear medium, 
reduces the cross section. This reduction is large in the energy region of Ey= 0.4-0.5GeV than the reduction in the 
energy region of -Ejy=l-3GeV as compared to the cross section calculated for the free case. However, in the case 
of antineutrino this reduction is more than in the case of neutrino as compared to the cross section calculated for 
the free case. In the incoherent charged current lepton production process accompanied by a pion, the reduction in 
the cross section due to the nuclear medium effects is more than the reduction due to the final state interaction of 
pions. While in the case of coherent pion production, the reduction due to final state interaction is quite large as 
compared to the reduction due to the nuclear medium effects. Furthermore, we find that the contribution of the cross 
section calculated in the case of coherent pion production process in the nuclear medium and final state interaction 
effects is around 6-7% to the total (incoherent+coherent) one pion production cross section in the energy region of 
QAGeV < E < 3GeV, and due to this we have not considered the form factor dependence, distribution, etc. in 
this case. 

(ii) We find that in the case of neutrino, the distribution calculated in the local Fermi gas model with RPA 
effects results in a large reduction in the peak region of Q^. For antineutrino induced process, there is a shift in the 
peak region which is towards low and the reduction is smaller than in the case of neutrino induced process. In the 
case of incoherent one pion production process, there is a reduction in the peak region of by taking medium effects 
into account which becomes smaller at large value of Q^. When pion absorption effects are also taken along with 
medium effects then there is a further reduction in the distribution. For antineutrino induced process, peak shifted 
towards low and the reduction is smaller than in the case of neutrino induced process. 

(iii) The dependence of the axial dipole mass Ma on the total scattering cross section for the neutrino induced 
charged current quasielastic process and for the incoherent charged current lepton production process in ^^O is 
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studied. We find that the scattering cross section increases when Ma is taken as 1.21GeV, while it decreases when 
Myi=1.05GeV as compared to the cross sections calculated with Ma— l.\Ge\ . In the case of antineutrino induced 
reaction the dependence on Ma becomes small. 

(iv) The dependence of the various parameterizations of the isovector form factors on the distribution in the 
case of charged current quasielastic lepton production process and N-A transition form factors in the case of charged 
current l7r"'"(l7r~) production process induced by neutrino (antineutrino) in ^^O has been studied. For the quasielastic 
process results have been presented with the parameterizations given by Budd et al. (5^1, Bradford et al. [s^, Bosted 
et al. [5^ and Alberico et al. 59]. Wc find that the use of various parametrization for the isovector form factor in the 
case of charged current quasielastic lepton production process results in a very small change in the peak region of . 

In the case of incoherent one pion production process these results are obtained by using the N-A transition form 
factor parameterizations given by Lalakulich et al. [T^I, Paschos et al. [7l| and Schreiner et al. [69]. We find that in 
the case of neutrino induced pion production process the differential cross section obtained by Paschos et al. [71] and 
Schreiner et al. [i^ is smaller in the region of low_Q^, which increases by small amount for Q^=0.2-0.4GeV^, than 
the cross section calculated with Lalakulich et al. 7^ parameterization. In the case of antineutrino induced process 
there is very small change in the distribution. 

(v) Nuclear medium effects play an important role in reducing the number of events obtained by integrating a over 
the atmospheric neutrino flux. In the case oiv -\-v induced quasielastic lepton production the reduction is around 25% 
when the events are calculated in the local Fermi gas model without the RPA effects and a total reduction of around 
45% when RPA effects are also taken into account in comparison to the lepton events calculated for the free case. In 
the case of inelastic lepton production the reduction from the events calculated without medium effects is around 40% 
in the case of incoherent pion production and around 85% in the case of coherent pion production, which results in 
a net reduction of around 50% when nuclear medium and final state interaction effects are taken into account. The 
lepton events from the inelastic process also contribute to the quasielastic events even when there are no pions in the 
final state because of its absorption while coming out of the nucleus. Such quasielastic like events contribute around 
12% to the total quasielastic lepton events. 

We find that for the total lepton events obtained by integrating the total cross section cross section(cr) over the 
atmospheric neutrino fiux results in a large reduction in the event rates when RPA effects are taken into account. 
Therefore, the results with nuclear medium effects for muon and electron events may be important in the analysis 
of neutrino oscillation experiments. Thus, we conclude that in the future neutrino oscillation experiments to be 
performed at SuperK on the atmospheric neutrino or the accelerator neutrinos like T2K and NOi^A, the study of 
nuclear medium effects in predicting the event rates would play a very important role. 
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